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Atmospheric Neutrino Physics with the MACRO detector 

M.Spurio a for the MACRO collaboration 

a Dipartimento di Fisica dell'Universita and INFN, 40126 Bologna, Italy 
e-mail: spurio@bo.infn.it 

We present the measurement of the the flux and angular distribution of atmospheric using the MACRO 
detector. Three different event topologies are detected in two different energy ranges. High energy neutrinos 
(E v ~ 80 GeV) via the identification of upward throughgoing muons. Lower energy neutrinos (E v ~ 4 GeV) 
via the upgoing stopping and partially contained downgoing muons (ID+UGS), or via the partially contained 
upgoing muons (IU). The measured flux is reduced with respect to the predictions. For the high energy sample, 
globally the flux reduction is 0.74 ± 0.054:stat+sys ± 0.12th; the reduction varies with the zenith angle. The 
ratio of measured to expected events is almost constant with the zenith angle for the low energy events, and is 
0.57±0.08 stt , t+s?/s ±0.14tw for the IU sample, and 0.71±0.09 st(1 i+ SJ/a ±0.17 t h eor for the (ID+UGS). All the data 
sets are consistent within a scenario of neutrino oscillations, with maximum mixing and Am 2 ~ 10 -3 -r- 10" 2 eV 2 . 



1. INTRODUCTION 

The interest on the atmospheric neutrinos has 
grown up in the last year, after the Neutrino '98 
Conference in Takayama, Japan. New, higher 
statistic data have been presented there by the 
Soudan 2 (I), MACRO g], and SuperKamiokande 
(SK) H collaborations. The measured flux of 
muons induced by atmospheric shows a reduc- 
tion with respect to the expectation, which de- 
pends on the neutrino energy and direction. For 
v e induced electrons there is no strong deviation 
from the prediction. The three experiments ex- 
plain the reduction in terms of neutrino os- 
cillations, with maximum mixing and Am 2 few 
times 1CP 3 eV 2 , confirming the early results. In 
fact, in the simplest scenario of two flavor oscilla- 
tions, the survival probability of a pure beam 
is: 



P{vy. —n/ tl ) = l — sin 2 28 sin 2 ( 



1.27Am 2 • L, 

El ' 



(1) 



Am 2 is the mass difference of the two neutrino 
mass states, 6 is the mixing angle, E v the neu- 
trino energy and L the path length from the pro- 
duction point to the detector. L can be estimated 
through the neutrino arrival direction O. For up- 
going neutrinos, as the zenith angle O changes, 
L ~ 2i?0 • cosQ (Rq is the Earth radius), while 
L is only few tens of kilometers for downgoing 



neutrinos. 

Atmospheric neutrinos are detected in the SK 
water Cherenkov detector via their interaction 
with water nuclei. Three different classes of 
events are defined (with increasing average en- 
ergy of the parent neutrino): fully contained 
events (FC), partially contained events (PC) and 
upward-going muons. Electron neutrinos are also 
identified in the FC sample. The ratio of muons 
to electrons normalized to the respective Monte 
Carlo predictions enhances the anomaly. 

The Soudan 2 results support the oscillation 
hypothesis by measuring atmospheric and v e 
interactions in the (roughly) same energy re- 
gion of SK. A different detection technique (drift 
chamber calorimeter) is used in this case. 

Here we present the MACRO results on the 
measurement of the atmospheric neutrino flux in 
the energy region from a few GeV up to a few 
TeV. In this case, a completely different exper- 
imental technique is used. The flux of is in- 
ferred from the measurements of upward through- 
going muons produced via charged current in- 
teractions in a large rock volume below the de- 
tector. The muon can travel up to the appara- 
tus, and if the residual energy is at least 1 GeV, 
it is detected as an upgoing muon. The average 
energy of the parent neutrino for these events is 
~ 80 GeV, one or two order of magnitude larger 
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than that of FC and PC Superkamiokande events. 
The muon can be detected in MACRO as an up- 
going stopping particle, if its residual energy is 
between 0.3 -f- 1 GeV. In this case, neutrino in- 
teraction happens few meters below the appara- 
tus, and the average parent neutrino energy is of 
the order of few GeV. In addition to the stop- 
ping muons, MACRO measures the flux of lower 
energy [E v ~ 4 GeV) neutrinos through the de- 
tection of (mainly) interactions inside the ap- 
paratus. 

2. MACRO AS v„ DETECTOR 
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Figure 1. Sketch of different event topologies 
induced by muon neutrino interactions in or 
around MACRO. The dark boxes represent the 
liquid scintillator hits (measurement of time and 
dE/dx), while the stars represent the hits on the 
tracking system. 



The MACRO detector is a large rectangular 
box (76.6 m x 12 m x 9.3 m) whose active de- 
tection elements are planes of limited streamer 
tubes for tracking and liquid scintillation counters 
for fast timing. It is located at the Gran Sasso 
Laboratory, with a minimum rock overburden of 



3150 hg/cm 2 . The lower half of the detector is 
filled with trays of crushed rock absorber alter- 
nating with streamer tube planes, while the upper 
part is open. The angular resolution for muons 
achieved by the streamer tube system is better 
than 1°. The time resolution of each liquid scin- 
tillation counter is about 0.5 ns. Fig. |l| displays 
the different kinds of measured neutrino events. 
The up throughgoing muons come from in- 
teractions in the rock below the detector. The 
muon crosses the whole detector and the flight 
direction is determined by time-of-fhght (t.o.f.) 
measurement. The flux of lower energy is stud- 
ied by the detection of interactions inside the 
apparatus; the partially contained upgoing events 
(IU) are tagged with t.o.f.. The partially con- 
tained downgoing events (ID) and upward going 
stopping muons (UGS) are identified via topo- 
logical constraints. Fig. || shows the distribution 
of the parent neutrino energy for the three event 
topologies detected by MACRO. The data pre- 
sented here come mainly from the running period 
with the full MACRO detector (which started ac- 
quisition in April 1994), up to February 1999, cor- 
responding to an effective live-time of 4.1 years. 




Figure 2. Expected distribution (for one year of 
data) of the parent v energy for the three event 
topologies in MACRO. 
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3. UP THROUGHGOING MUONS 
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Figure 3. 1/(3 distribution (data with the full 
detector); the upgoing muons are centered at 
1/(3 = — 1. The shaded distribution is for the 
subset of events with three scintillator layers. 



At the depth of the Gran Sasso Laboratory, we 
expect ~ 5 x 10 4 downgoing atmospheric muons 
for each neutrino-induced /i (for which we ex- 
pect an up/down symmetry). For this reason, 
the identification of neutrino-induced /i's relies 
on the measurement of the direction that muons 
travel through MACRO. For each detected muon 
the experimental parameter 1/(3 — c '( Tl ~ T2 } is 
evaluated. Tj. (T2) is the time measured in the 
lower (higher) scintillation counter and D the 
path length between the counters. The downgo- 
ing muons are expected in the 1/(3 region near 
+1, while upgoing muons at 1/(3 near -1. To 
remove accidental background events, the posi- 
tion along the scintillator counter must agree 
within ±70 cm with the position indicated by 
the streamer tube track. Fig. || shows the dis- 
tribution of 1/(3 for muons collected from April 
1994. There are 561 upgoing muons in the range 
-1.25 < 1/(3 < -0.75. We combine these data 
with additional 81 events collected before 1994, 



for a total of 642 upgoing events. In the to- 
tal data set there are: 12.5 ± 6 estimated back- 
ground events due to misidentification of downgo- 
ing muons; 10.5 ± 4 background events from up- 
going charged particles produced by downgoing 
muons in the rock near MACRO 12 ± 4 inter- 
nal events from interactions of neutrinos in the 
very bottom layer of MACRO scintillator. Re- 
moving the background and the internal events, 
the number of upward throughgoing muons is 607. 

3.1. Monte Carlo expectation 

For the simulation, we have used the Bartol 
neutrino flux [g and the DIS parton distribution 
set for the neutrino cross-sections. The prop- 
agation of muons to the detector has been done 
using the energy loss calculation for standard 
rock. The total uncertainty on the expected flux 
of muons adding in quadrature the errors from 
neutrino flux, cross-section and muon propaga- 
tion is ±17%. This theoretical error in the predic- 
tion is mainly a scale factor that doesn't change 
the shape of the angular distribution. The num- 
ber of expected events is 824.6, giving a ratio of 
the observed number of events to the expecta- 
tion of 0.74 ± 0.031 stat ± 0.044 sys ± 0.12 theo . Fig. 
|] shows the zenith angle distribution of the mea- 
sured flux of up throughgoing muons (all MACRO 
data), compared to Monte Carlo expectation. 

3.2. Interpretation of the result 

We interpreted the reduction on the detected 
number of events and the deformation of the 
zenith angle distribution as a consequence of 
disappearance. In the scenario described by eq. 
1, relatively fewer events are expected near the 
vertical (cosO = —1) than near the horizontal 
(cosO = 0), due to the longer path length of 
neutrinos from production to observation. The 
shape of the angular distribution has been tested 
with the hypothesis of no-oscillation. We found a 
x 2 /d.o.f = 22.9/8, or a probability of 0.35%. To 
test the — > v T oscillation hypothesis, we evalu- 
ate the independent probability for obtaining the 
number of observed events and the angular dis- 
tribution for various oscillation parameters. The 
maximum of the \ 2 probability in the physical 
region of the oscillation parameters is 36.6%, cor- 
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Figure 4. Measured flux of upward throughgoing 
muon vs. the cosine of zenith angle 9. The 17% 
uncertainty (shadow) is almost a constant correc- 
tive factor to the central value, being the error on 
the shape almost negligible. The lower line shows 
the prediction assuming two-flavors neutrino os- 
cillations (see text). 



responding to Am 2 = 2.5 x 10~ 3 eV 2 and maxi- 
mum mixing. Fig. || shows the confidence regions 
at the 90% and 99% C.L. in the parameter space 
(sin 2 29, Am 2 ) for — > v T oscillations, based on 
application of the Monte Carlo prescription of || . 
The effect of oscillations (using the best fit point 
values) is also shown in Fig. ||. The angular dis- 
tribution of events vs. cosQ strongly disfavours 
the hypothesis of no-oscillations; however, there 
is a structure (near cosQ ~ —0.65) which is unex- 
pected also in the case of oscillations. The excess 
of events in this bin is still consistent with a sta- 
tistical fluctuation, but it could be a hint for a 
more complex scenario. 

3.3. Checks on systematics 

Possible systematic effects have been studied 
and have been shown to be too small to ex- 
plain the observed anomalous shape in the zenith 
angle distribution. The detector acceptance is 
well understood (from a large sample of down- 




stage 

Figure 5. Allowed contours (and best fit point) at 
90% and 99% C.L. (darker region) assuming — » 
v r oscillations. The 90% C.L. contour obtained 
from the lower energy neutrino events (IU and 
ID+UGS) is also shown. 



going muons), in particular near the vertical di- 
rection where the largest deviation compared to 
the Monte Carlo expectation without oscillations 
is observed. 

An important check has been performed using 
the system designed for the detection of bursts of 
neutrinos from stellar collapse. This system has a 
completely separated acquisition system and red- 
out electronics. An analysis has been done using 
the TDC's of such system and the events having 
three scintillation counters. With three counters 
is possible to achieve a good rejection of downgo- 
ing muons using only the information from this 
TDC system and without tracking. The angular 
distribution of the up throughgoing event sample 
identified in this way (presented in Fig. ^]) con- 
firms the main feature of Fig. || In the figure 
some events collected when the stellar collapse 
system was on and the main system off are in- 
cluded. During the period of common activity, 
only one event has escaped from the detection of 
the standard analysis, while we expect to miss 
0.6 events due to dead time of the the streamer 
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tube/liquid scintillation systems used by the stan- 
dard analysis. 
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Figure 6. Zenith angle distribution of the up- 
going muons data (points) compared with the 
Monte Carlo predictions without oscillations for 
the events measured with the system dedicated 
for the search of stellar collapses (~ 1.7 y live 
time) 

4. EVENTS FROM LOWER ENERGY 

4.1. Partially contained upgoing /i (IU) 

The detection methods for the IU are similar to 
that of throughgoing /i's (time-of-flight measure- 
ment, plus a track reconstructed in the streamer 
tubes system), and on topological criteria for the 
identification of a interaction vertex inside the ap- 
paratus. To reject fake semi-contained events en- 
tering from a detector crack, the extrapolation of 
the track in the lower part of the detector must 
cross and not fire at least three streamer tube 
planes and one scintillation counter. The above 
conditions, tuned on the Monte Carlo simulated 
events, account for detector inefficiencies and re- 
duce the contribution from upward throughgoing 
muons which appear like semi-contained to less 
than ~ 1%. We evaluated that 5 events are due 
to an uncorrelated background. After the back- 
ground subtraction, 116 events are classified as 
IU events. 



4.2. Partially contained downgoing and up 
stopping events (ID) 

The ID+UGS events cross only one liquid scin- 
tillator layer and are identified by means of topo- 
logical criteria. The lack of timing information 
prevents to distinguish between the two sub sam- 
ples. From MC (sec. 4.3.), an almost equal num- 
ber of UGS and ID events are expected. A soft- 
ware selection rejects throughgoing events, and 
searches for a contained track crossing the bot- 
tom layer of the scintillation counters. The se- 
lection conditions for the event vertex (or fi stop 
point) in the detector are symmetrical to those 
for the IU search. 879 events are accepted by the 
selection. Some of them are wrongly tracked or 
bending atmospheric muons which entered from a 
detector crack. To reject fake events, a visual scan 
was performed. Two physicists scanned twice the 
real data randomly merged with the simulated 
events. At the end of the scan procedure, 200 
real events are accepted as ID or UGS, with 95% 
of the Monte Carlo simulated events. The main 
background source are upward going charged pi- 
ons induced by interactions of atmospheric muons 
in the rock around the detector. 7.2 ± 2.3 back- 
ground events have been evaluated using a full 
simulation which is based on our measurement 

i- 

4.3. Monte Carlo 

Also the expected low energy event rates have 
been evaluated with a full Monte Carlo simula- 
tion. The interaction of atmospheric v e and j^] 
were simulated in a large volume (with 175 kton 
total mass), including the experimental hall and 
the detector. Because of the lower v energies, the 
cross sections from [Q are used, which include the 
contribution of the exclusive channels of quasi- 
elastic scattering and single-7r production. The 
total theoretical uncertainty on v flux and cross 
section at these energies is ~ 25%. The detec- 
tor response has been simulated using a GEANT 
based program, and events are processed in the 
same analysis chain as the real data. The parame- 
ters of the streamer tube and scintillator systems 
have been chosen to reproduce the real average 
efficiencies. A 10% systematic error is evaluated 
from the simulation of detector response, data 
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taking conditions, analysis algorithm efficiency, 
mass and acceptance of the detector. As shown 
in Fig. ^, the energy spectra of parent neutrinos 
for the IU and ID+UGS events are similar, with 
equal average energy. 87% of IU and (ID+UGS) 
detected events are induced by z^-CC interac- 
tions. The remaining 13%, from z^ e -CC and NC 
interactions (with different percentage for the two 
data sets). 

4.4. Results from the low energy events 

Fig. ^ shows the zenith angle distribution of 
the IU and U GS + ID data samples, with the 
Monte Carlo predictions. The data are within 
errors consistent with a constant deficit in all 
bins with respect to the Monte Carlo expecta- 
tions. The ratios of the number of observed to ex- 
pected events are Rid+ugs = (^^ l )id+ugs = 
0.71 ± 0.05 stat ± 0.07 SJ/st ± 0.17 theor and R w = 
0.57 ± 0.05 stat ± 0.06 syst ± O.U theor . 

If the event deficit is due to an overall theoreti- 
cal overestimate of the neutrino flux and/or cross 
sections, it is expected 1Z = Rw / Rid+ugs = 1. 
The theoretical and systematic errors are largely 
reduced (to 4% and 5%, respectively) if the ratio 
of ratios is considered. The partial uncertainty 
cancelation comes from the fact that the energy 
spectra of the two topologies are alike, and be- 
cause of the symmetry in the detector acceptance. 
We measured 1Z = 0.80 ± 0.09 stat ; the statisti- 
cal error is the dominant one in this quantity. 
The probability to obtain a ratio so different from 
the expected one is 5% (taking into account the 
non-gaussian shape of the uncertainty) , no matter 
which neutrino flux and neutrino cross sections 
are used for the predictions. 

The alternative hypothesis for such a reduction 
is fn — > v T oscillations with maximum mixing 
and Am 2 ~ 10 -3 10~ 2 eV 2 . In this Am 2 in- 
terval the number of muons induced by the in- 
teraction of upgoing neutrinos (i.e. IU and UGS 
events, for which L ~ 13000 km , E v ~ 4 GeV) 
is reduced by a factor of two, while almost no re- 
duction is expected for ID events (L= few km). 
The expected (ID+UGS) event rate is 3/4 of the 
no-oscillations expectation (neglecting the contri- 
bution of v e and NC interactions). For larger 
Am 2 , also the ID events are reduced, so both the 



(ID+UGS) and IU event rates are 1/2 of the no- 
oscillations expectation. For smaller Am 2 , the 
bins of the zenith distributions (Fig. 0) are dif- 
ferently affected by the oscillations. As the low- 
energy events are not particularly sensitive to the 
Am 2 interval ~ 10~ 3 -j- 10~ 2 eU 2 as a test point 
in Fig. ^ we used the best fit values from the high 
energy sample. 

Finally, the — > v T hypothesis simply ex- 
plain the different reduction on the total num- 
ber of events observed in the (ID+UGS) and IU 
data sets. In case of v„ disappearance, for the 
(ID+UGS) the reduction with respect to the case 
of no-oscillations is 0.76, to be compared with the 
measured value of Rid+ugs = 0.71. For the IU 
events, it is 0.57, to be compared with Rw — 0.57 

We estimated the most likely values on a 
(sin 2 28, Am 2 ) grid using a \ 2 comparison of data 
and Monte Carlo, based on the prescription of ||. 
The data were binned in 4 zenith angle bins for 
the IU events, 4 zenith angles for the ID+UGS 
events, the ratio id+iigs anc ^ ^ e overau nor ~ 
malization. The maximum of the x 2 probability 
(97%) occurs at sin 2 26 = 1.0 (inside the physi- 
cal region); this value of the \ 2 probability is al- 
most constant in the interval Am 2 = 10~ 3 + 2. x 
10~ 2 eU 2 . Fig. H shows the contour 90% confi- 
dence level for the low energy events; the allowed 
region is consistent with that obtained using the 
higher energy sample of neutrino-induced upward 
throughgoing muons. 

5. CONCLUSIONS 

MACRO measures three different data sample 
of events induced by atmospheric neutrinos. All 
the data sets show a deficit of the measured num- 
ber of events with respect to the predictions based 
on the Bartol flux and the absence of neutrino os- 
cillations. 

For the upward throughgoing sample, apart 
from the reduction in number, the shape of the 
angular distribution is modified as expected in 
case of disappearance. There is however an 
unexpected structure near cosQ ~ —0.65), which 
is still consistent statistical fluctuation in 

the scenario of neutrino oscillations. 

For the low energy neutrinos, two samples are 



7 



c 120 
A 00 

0) 

40 



20 




UGS+ID 



-0.6 -0.4 



-0.2 

— IcosGI 



120 
100 



2 80 

| 60 
CD 

40 
20 




I J 



1 -O.i 



-0.6 -0.4 



-0.2 
COS© 



Figure 7. Cosine of the zenith angle (6) dis- 
tribution for (ID+UGS) and IU events. The 
background-corrected data points (black points 
with error bars) are compared with the Monte 
Carlo expectation assuming no oscillation (full 
line) and two-flavour oscillation (dashed line) us- 
ing maximum mixing and Am 2 = 2.5 x 1CP 3 eV 2 . 



measured. The IU events are induced by up- 
going neutrinos, while the ID+UGS events are 
50% from neutrino from above (ID). The mea- 
sured number of events is below the expecta- 
tions, but the reduction is different for the two 
data set. Because the parent neutrinos for the IU 
and the ID+UGS events have a quite similar en- 
ergy spectrum, an overall reduction of the num- 
ber of neutrino-induced muons has a low prob- 
ability (5%) to explain the two observed differ- 
ent deficit. This effect is explained with larger 
probability (~ 97%) by the hypothesis of muon 
neutrino oscillations with maximum mixing and 
Am 2 = 10" 3 v2x 1CT 2 eU 2 . 

In Table 1 is presented a summary of the 
number of detected events, the number of the 
expected ones in case of no-oscillation and for 
—> v t oscillations with maximum mixing and 
Am 2 = 2.5 x 10~ 3 eU 2 . The up through-going 
sample has the higher sensitivity to the oscilla- 



tion parameters. For this sample, the detection 
technique and the average energy of the parent 
neutrino are completely different from the SK 
and Soudan 2 ones. However, the allowed region 
in the oscillation parameter space largely overlap 
that of SK and Soudan 2, and also the best fit 
parameters are similar to the SK one. 





Data 


MC (Bartol flux) 






No osci With osci 




- Bckg 


± s ta+ sys 0.0025eV 2 


Up Through 


607 


824 ± 56 585 


IU 


116 


202 + 22 115 


ID+UGS 


193 


273 ± 30 209 



Table 1 

Summary for the MACRO measurements of the 
atmospheric v flux. For the three event topolo- 
gies, the number of detected events (first column) 
is compared with the expectation in case of no- 
oscillations, and for za, — * v T oscillations with 
maximum mixing and Am 2 = 2.5 x 10 -3 eV 2 
(best fit point for the high energy set). 
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